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ABSTRACT 

Characterizing the conversion factor between CO emission and column density of molecular hydro- 
gen, Xco, is crucial in studying the gaseous content of galaxies, its evolution, and relation to star 
formation. In most cases the conversion factor is assumed to be close to that of giant molecular clouds 
(GMCs) in the Milky Way, except possibly for mergers and star-bursting galaxies. However, there 
are physical grounds to expect that it should also depend on the gas metallicity, surface density, and 
strength of the interstellar radiation field. The Xqq factor may also depend on the scale on which 
CO emission is averaged due to effects of limited resolution. We study the dependence of Xco on 
gas properties and averaging scale using a model that is based on a combination of results of sub-pc 
scale magneto-hydrodynamic simulations and on the gas distribution from self-consistent cosmological 
simulations of galaxy formation. Our model predicts Xco ~ 2 — 4 x 10^° K^^ cm^^ km~^ s, consis- 
tent with the Galactic value, for interstellar medium conditions typical for the Milky Way. For such 
conditions the predicted Xco varies by only a factor of two for gas surfaced densities in the range 
^ 50 — 500 Mq pc^^ . However, the model also predicts that more generally on the scale of GMCs, 
Xco is a strong function of metallicity, and depends on the column density and the interstellar UV 
flux. We show explicitly that neglecting these dependencies in observational estimates can strongly 
bias the inferred distribution of H2 column densities of molecular clouds to have a narrower and offset 
range compared to the true distribution. We find that when averaged on ~ kpc scales the X-factor 
depends only weakly on radiation field and column density, but is still a strong function of metallicity. 
The predicted metallicity dependence can be approximated as Xco oc Z^^ with 7 sa 0.5 — 0.8. 
Subject headings: galaxies: evolution - galaxies: ISM - ISM: molecules - methods: numerical 



1. INTRODUCTION 

Molecular hydrogen (H2), the major constituent of cold 
clouds in the interstellar medium (ISM) , is playing a ma- 
jor role in shaping the visible Universe around us. Cool- 
ing by H2 aided the formation of the first stars (e.g., 
lAbel et al.ll2002l : iBromm et alll2002[ ) and star formation 
in nearby galaxies takes place in molecular and giant 
molecular cloudf0 (GMCs). Furthermore, the proper- 
ties of molecular gas in galaxies, e.g., its distribution, 
formation and destruction are intricately linked to many 
other phenomena occurring in galaxies, such as the for- 
mation of stars, various stellar feedback processes, and 
the physics of dust grains. 

Unfortunately, detecting H2 in emission is difficult be- 
cause the lowest excited levels in the rotational ladder are 
hardly populated at the low temperatures (~ 10 K) of 
molecular clouds. Moreover, the H2 molecule, lacking a 
permanent electric dipole moment, radiates via the much 
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slowe r quadrupole transit ion (e.g., IShull fc BeckwithI 
1982; Stabler fc PalIall2005D . Therefore, a practical solu- 
tion has been to measure the emission from some other 
molecule, that is assumed to trace H2, and convert its 
line intensity into an H2 column density. The conversion 
factor is called the X-factor for the particular emission 
line. 

Cooling lines from the carbon monoxide (CO) iso- 
tope ^^CO are often used as H2 tracers, because 
^^CO is the second most abundant molecule in molec- 
ular clouds and can be observ ed even in extragalac- 
tic o bj ects at high redshifts ([Brown fc Vanden Bout] 
1991"; 'Solomon & Vandcn Bout' '2005'; 'Tacconi et aU 
2006,. 2008: Daddi et al. 2010a,b- Tac_C(mi et al, 201fl 
Genzel et al.l 120101: lEmonts et al] l20linGenzel et alJ 
201 Ih . In this study we will focus on the ^^CO J 



1—^0 transition (and on the corresponding X-factor 
Xco) I i-fi-i the transition from the first excited rota- 
tional level to the ground state, that is widely used in 
the literature, particularly in studies of nearby galax- 
ies or molecular clouds in the Milky Wa v (I Wilson et aLl 



1970; Scoville & Sanders 1987; Solom on et al.l 11987 



Young & S covilla 1I99I1: lY'oung et al.l [1995 1: 'Rega n et al 
|001; Hei fer et al.l 120031 : iKuno et al.l 12007: Blitz et al 
120071) . 

Xco has been measured for molecular clouds in the 
Milky Way and in a few other nearby galaxies with a va- 
riety of techniques. These methods compare the directly 
measured CO luminosity (or intensity) with an indepen- 
dently determined H2 mass (or column density). The 
most common approaches: (i) estimates of virial masses 
from line- width and cloud size (e.g., iSolomon et al.l 
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I1987D . (ii) measurement of H2 masses using a different 
tracer molecule (e.g., ^■^CO), or a higher order t ransi- 
tion, with its ow n X-factor (e.s.. IDickm;ml [19751 . \WM 
iHever et al.l 120091 ) . (iii) measurement of proton number 
density from gam ma rays produced in cosmic ray - pro- 
ton c o llision (e.K..lBloemen et al.lll986l:lStrong fc Mattoxl 
[1991 IHunter et al.l 119971: lAbdo et al.ll20Tfl[ i. Cw) mea- 
surement of to tal hydrogen column density using infr ared 
emission fe.g.. lDame et al.ll2001l : [Praine et al.ll2007| ). (v) 
estimates of visual extinction. Ay, and thus total gas 
density from IR sta r counts (l^yol^l 119231) or the near- 
infrared color excess ()Lada et al.lll994l ). The latter three 
approaches have to be supplemented with maps of atomic 
hydrogen unless most of the hydrogen is molecular. 

The consensus is that for GMCs in the Milky Way 
all these different methods give similar values of Xqo ~ 
1.5 — 4 X 10^" cm~^ km~^ s. The question why the 
X-factor is remarkably constant in molecular clouds in 
the Milky Way is not yet entirely settled, but it has been 
suggested to be a consequence of the narrow range of H2 
column densities and temperatures of Milky Way clouds 
(IShettv et al.l[20Tlal) . 

On the other hand, there is increasing evidence 
that the X-factor may be different in othe r galaxies 
(e.g., iBolatto et all 120081 : iLerov et al.l 120 lit ). Specif- 
ically, there has been a long debate as to whether 
the X-fac tor i s larger in low metallicity galaxies (e.g., 
Wilson ,1991 lArimoto et"ai.. 1996: Boselli ct al. 2002; 
Israeli 120051 but cf. iBlitz et al.l I2007D . and possi- 



bly s maller in the central regions of disk galaxies 
(e.g., iMalonev fc Blackl 119881 : lOka et al.lll998[ ). in local 
starbursts and (ultra-) luminous in frared galaxie s (e.g.. 
Wild et al.l 1992; Shier et al.l 119941: i Maucrsbcrg er et alJ 



19961: iSolo mon ct al. 1997; Downes & SolomonnTgM 



Brvant fc Scovill e 1999: Meier et al... 20101) or insubmil 



limeter galaxies at higher redshifts (e.g., iTacconi et al.l 
|2008| ). In a simple toy model in which the CO emis- 
sion stems from an ensemble of optically thick. 



alized clumps, the X-factor should scale as 



''H 



IT 



([Dickman et al.l [19861) . However, one arrives at a more 
complex scaling of the X-factor if the complex geome- 
try of the supersonically turbulent molecular gas and ef- 
fects of radiative t ransfer are taken into account (e.g., 
IShettv et al.ll2011b[ ). 

A further important aspect is that many extragalactic 
surveys do not resolve individual GMCs, but measure CO 
emission on ~ kpc or even larger scale. A study of the 
properties of the X-factor on such large scales, however, 
requires galaxy models with reasonably realistic density 
distributions. The large dynamical scale necessary to re- 
solve individual CO emitting gas clumps ('-^ 0.1 — 1 pc) in 
a self-consistent, preferentially cosmological, simulation 
poses a serious challenge. 

In this paper, we study the behavior of the X-factor as 
a function of gas metallicity, UV radiation field, and as a 
function of scale. We analyze the impact of the X-factor 
on star formation laws in a follow-up paper (Feldmann 
et al, in prep). In section [2] we present our novel nu- 
merical approach that aims at circumventing the prob- 
lem of the large dynamical scale. The basic idea is that 
we combine cosmological simulations, in post-processing, 
with a model that is calibrated on results of small-scale 
magneto-hydrodynamic (MHD) simulations of the ISM. 



This approach allows us to study CO emission at GMC 
scales and above with an effective resolution of ~ 0.1 pc 
in a cosmological volume. In section 13.11 we discuss the 
scaling of the X-factor with II2 column density on GMC 
scales. Next, in section l373l we study the dependence of 
Xqo on the local environment (metallicity and UV ra- 
diation field) . Subsequently, in section 13.41 we discuss 
the observational evidence of a constant surface density 
of molecular clouds. Finally, in section |4l we address 
the effect of spatial averaging of the X-factor. We then 
summarize our results and conclude. 

2. METHODS 

The approach we adopt in this study is to complement 
large scale simulations, in which the global structure of 
galaxies and their ISM is modeled self-consistently in a 
cosmological context, with an adequate subgrid model 
for CO emission. 

One possibility is to c ompute the X-factor based 
on PDR models (e .g. iTielens fc HollenbachI [l985l: 
van Dishoeck fc Blackl Il988t [Sternberg fc Dalgarnol 



1995^ 'Kaufman et^ 120061 : iWolfire et al.l I2OIOI : 



Krumholz et al. 201l[). Although the PDR approach is 
useful in constructing simple models to make predic- 
tions about the atomic-to-molecular transition or the 
X-factor, one may wonder how sensitive such predictions 
are to the specific assumptions of this method. It is 
thus useful to explore alternative models based on a 
different set of assumptions. To this end, we use a model 
based on tailored small-scale simulations that are able 
to res olve the turbulent structu re of the ISM on ^ pc 
scales (jGlover fc Mac Lowl[20TT[) . 

Doing so has several key advantages. First, these small 
scale simulations rely on well-defined physically under- 
stood processes and their incorporation into cosmolog- 
ical simulations fits well within an 'ab initio' approach 
of understanding the evolution of galaxies and the ISM. 
Second, the small scale simulations follow the chemical 
evolution and the dynamics of the ISM in a presumably 
more realistic way than ad hoc sub-grid models. Fi- 
nally, large scale simulations can provide the appropriate 
boundary conditions for the small scale simulations and 
thus, as long as one stays within the regime explored 
by the small scale simulations, this approach is poten- 
tially not plagued by free parameters that have to be 
tuned. Thus, if the model turns out to disagree with ob- 
servations, we cannot simply adjust the parameters of the 
model. Instead, this would signal that important physics 
is missing which needs to be included in the small scale 
simulations. 

However, there are a few technical challenges to this 
op timistic perspective . First, since the ISM simulations 
of [Clover fc Ma c Low (2011) assume a Milky- Way like 
UV radiation field (Draine 1978t) . we have to model the 
UV dependence of the CO abundance. We do this using 
a simple but approximate assumption that CO is in pho- 
todissociation equilibrium. Second, we do not perform 
line radiative transfer (RT). Instead, we use the photon 
escape probability formalism tha t has been shown to pro- 
vide a good fit to observations ([Pineda et al.l 120081 ) and 
compares reasonably wel l with line RT calcu lations, but 
may differ in some detail ([Shettv et al.|[2011al) . Third, we 
can neither model the CO line-width nor the excitation 
temperature self-consistently at our numerical resolution. 
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We address the first point by considering two extreme 
cases of the scaling of the line width with column density 
finding some quantitative, but not qualitative, changes in 
our results. We assume that the excitation temperature 
is constant, which is roughly correct for many GMCs 
in the MW and moderately star fo rming galaxies (e.g., 
iLerov et"aLll2009HHeYer et al.ll2009[ ). but this assumption 
may break in galaxies undergoing massive starbursts or 
in high redshift mergers (see, e.g.. iNaravanan et al.ll20l"ll 
and references therein). Hence, while our approach al- 
lows a precise modeling of CO and H2 abundances, it 
does rely on assumptions regarding the gas temperature 
and the scaling of the velocity dispersion to infer CO 
intensities. This makes it an orthogonal ansatz, with 
its own strength s and weaknesses, compared with other 
appro aches (e.g.. lShettv et al.ll2011al |bl: INaravanan et al.l 
l2011f ). Many of these limitations may be circumvented 
in the future when a larger set of small-scale ISM simu- 
lations is available. 

A main limitation of our approach is that it relies on 
the assumption that the small scale simulations are a 
reasonable approximation to the real ISM, even though 
several physical mechanisms, such as star formation 
and stellar feedback, have not been included. Fur- 
thermore, possible feedback from small to large scales 
is not accounted for in our approach. Modeling the 
ISM below the scales of GMCs is a very challenging 
problem due to the variety of physical processes and 
feedback mechanisms that are operating such as stel- 
lar winds, HII gas pressure, protostellar jets or radia- 
tion pr essure (e.g.. iMatznen l2002t IMac Low fc KlessenI 
20041: iLTfc Nakamural 120061: iKrumholz et al.l I2Q06I: 



2009t iKru mholz & Matzner '2009: 
Murray et al. 2010 ; Diiajnc 2011 ; 



Nakamura fc Lil 120071: iBaneriee et al.| 120071: iBrunt et"aLl 



Wang ct al.' '2010^^ 

Lopez ct al. 2011; 

Murray et all I2OIII : fColdbaum et aP l201ll ). On su- 
per CMC (> 100 pc) scales feedback from su- 
pernovae and (potentially) active galacti c nuclei are 
likely important fe edback mechan i sms (e.g..lSilk fc ReesI 
1998; Mac Low fc KlessenI l2004t iFeruglio et al.l 120101 : 
Murray et al.il201lt ). 

A reasonable scale for separating sub- and super-grid 
physics lies therefore in the range ~ 10 — 100 pc and the 
resolution of our simulations (formal grid resolution is 
^ 60 pc) is chosen to fall in this window. The resolu- 
tion roughly matches t he box-sizes (20 pc) of the ISM 
simulation presented in iGlover fc Mac Lowl f2011) that 
we employ to extend the resolution of our simulations to 
smaller scales. 

Note that we use Zq = 0.02 (12 + logio(0/i?) = 
8.92) throughout the paper, the metallicity of the so- 
lar neighborhood, which is somewhat larger than the 
metallicity of the Sun according to recent estimat es 
(jAUende Prieto et al.l[2Mt IAsplund"enil[200l [20091 ). 

We present our X-factor model in ^2.1[ We give a short 
description of the suite of simulations that we use in this 
paper in In ii2.3l we discuss how we add the X-factor 
model as a post-processing step to our simulations. 

2.1. Modeling of the X-factor 

In this section we will discuss a simple model for the X- 
factor that is calibrated using small scale (~ few tens of 
parsecs) MHD simulations. The aim is to predict Xqo on 
these scales as a function of metallicity Z, mean density 



hh, and the UV radiation field Umw- The latter is given 
in units of the local int erstellar radiation field (|Drainel 
flgTllMathis et al.lfT98l . i.e. 

t^MW = •^ioooa/'^mw, 

where ^^loooA the mean intensity at 1000 A and Jmw — 
10^ photons cm~^ s^^ sr^^ eV^^. 

The X-factor is defined as the ratio between the molec- 
ular hydrogen column density iVjja and the velocity in- 
tegrated intensity of the J = 1 — > transition of carbon 
monoxide along a line of sight: 



X. 



CO 



Wco 



(1) 



A canonical value of Xco of molecular clouds in the 
Milky Way is Xqomw = 2 x lO^" cm'^ K"! km"! s, 
which we will refer to as the galactic X-factor. 

Computing Wco in a simulation is non-trivial unless 
the gas is optically thin to its own line radiation. We 
use an approximate treatment based on the escape prob- 
ability formalism to acco unt, in a crude fashiorQ, fo r line 
radiative transfer effects ([Glover fc Mac Low|[20lT[ ) 



T^co = TbAv 



2P{T)dT. 



(2) 



Here, Au is the width of the CO line, rio is the optical 
depth of the ^^CO J = 1 -» transit ion, ^(r) is the 
photon escape probability ([Tielensll2005D 



[1- 



,-2.34t 



]/(4.68r) if r < 7 



l/(4T[ln(r/V^)]i/2) ifr>7. 



(3) 



and Tb is the brightness temperature of the line, which 
we compute based on the gas temperature T and the 
CMB temperature Tcmb = 2.725(1 -I- z) as 



1 



1 



35.5/T _ I g5.5/TcMB _ 1 



(4) 



Despite the simplicity of the method, its predictions com- 
pare reasonably well with more elaborat e line radiative 
transfer calculations (iShetty et al.ll2011al) and it provides 
a good fit to observations ([Pineda et al.ll2008D . 



In local thermal equilibrium (LTE) the optical depth 
of the J = 1 ^ line is given by ([TielensI [200l 
IStahler fc Pallal[2005[ ) 



no ~ 1.4x10 



-16 



-5.5/T 



10 K 



T 



3 km s 



CO 



cm 

Given estimates of the CO and H2 column densities, and 
with some assumptions about the temperature and line 
width (see !j2.3l) . we can compute the X-factor using ([T]- 

U- , , 

Table 2 of IGlover fc Mac Lowl ([20111 ) contains the 
mass- weighted mean abundance of CO and H2 at the end 
of each of their simulations. The mass- weighted mean 
abundance Xs of species s (either CO or H2) in a box 
with total mass M = ^ p/XV is defined as 



Xs 



M 



E 



ns{i) 



p{i)XV(i) 



Us 



as 



71h(?) 

^ The approach assumes strictly local photon absorption and 
plan-parallel geometry. 
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Here, the index i runs over all grid cells in the box, while 
riH and Hs denote the local number densities of hydro- 
gen nuclei and molecules of species s, respectively. The 
convenient choice = 2 (while aco = 1) implies that 
Xfjs spans the range — 1 . The equivalent formulation on 
the right hand side gives the mass- weighted mean abun- 
dance in terms of the mean number density of molecules 
of species s, n^, and the mean number density of hydro- 
gen nuclei nn- 

We compute the average surface density Ng of species 
s along a line-of-sight of length L as 

N, = ^uhL. (6) 

The visual extinction along a line-of -sight is propor- 
tional to the column density of dust grains. Since the 
gas-to-dust ratio scales with the metallicity of the gas 
at least down to metallicities of ~ 1/10-th solar (e.g., 
iLerov et"aLll2011[ ). we compute the mean (averaged over 
lines-of-sight) visual extinction as 

Ay ^iz^ fiH L, (7) 
^0 

with 7z = 5.348 X 10^22 cm^ Z'^. 

Most of the simulations of iGlover fc Mac Low! (|2011f ) 
assume Umw = 1- However, the CO fraction in molecu- 
lar clouds should depend to some degree on the strength 
of the incident UV radiati on (e.g.. Ivan Dishoeck fc Blac5 
119881 iWolfire et al]l2010[ ). In highly star forming en- 
vironments, such as gas-rich galax y mergers or young 
galaxies at high redshifts (e.g., iStacev et al.l 119911 : 
llvison et al.|[2010l ). or in low dust galaxies such as low 
metallicity local dwarfs (e.g, iSmith fc Hancoclj l2009f ). 
the radiation field is expected to be much larger than 
that incident onto a typical GMC in the Milky Way. 
In order to proceed we will make a number of simpli- 
fying assumptions based on the simulation results of 
IGlover fc Mac Low! ()201lD (see their Fig. 1 fc 2). First, 
we assume that xu^ depends on uh, L and Z only 
through the combination uhZ , i.e. the number density 
of dust grains available to form H2. Specifically, we will 
assume 

xn^ {nn, Z, L, Umw) = {Zuh, C/mw)- (8) 

Several prescriptions exist in the literature to com- 
pute them^aiijn^lB^^ 

Il988t IKrumholz et"alll2008t IGnedin fc Kravtsov! 1201 iD . 
An alternative is to use the H2 abundance that i s esti- 
mated in the simulations of IGlover fc Mac Lowl ()2011[ ) 
and ignore the explicit radiation field d ependence which 
is ex pected to be only logarithmic (jKrumholz et al.l 
120081) . This latter approach is not as crude as it sounds, 
because the gas will be predominantly molecular when- 
ever there is significant CO emission and hence self- 
shielded from UV radiation. 

Second, we assume that xco depends mainly on the 
mean extinction and the radiation field: 

xcoiZ,nH,L,UMw) ^ xco{AviZnHL),Uuw, Z), (9) 

yet we have also included an explicit dependence of the 
CO fraction on metallicity, which is required because 
the CO abundance can never be greater than the abun- 
dance of C nuclei xc in the gas. We model this by 



enforcing an upper limit xc = 1.41 x 10^^ Z/Zq for 
xco{Av,Umw, Z). The crucial point is that the sim- 
ulations presented in IGlover fc Mac Low! (|20lH ) provide 
xco{Av,Umw = i,Z). 

In order to model the dependence of xcoiAy, Umw, Z) 
on the radiation field we assume equilibrium between 
the rates of photodissociation of CO and its creation by 
chemical reactions, i.e. 

a^co C^MW '5'dust(^v) '5'h2(^H2) "5*00(^^00) = nn'^TlijXi 

(10) 

where Sdust, Syi^, Sco are the shielding factor^] due to 
dust, H2 and CO, respectively, and the r.h.s. is a sum 
over all relevant chemical reactions that produce CO. 

We determine xco(^V7C^mw 7^ ^,Z) in the following 
way. Consider a box of size L, mean density fiH, metal- 
licity Z, and visual extinction Ay oc ZuhL exposed to a 
radiation field Umw 7^ 1- Consider now a second box of 
the same size L' — L, same metallicity Z' = but with 
possibly different mean density and thus visual ex- 
tinction Ay that is exposed to a MW-like radiation field 
t^MW = 1- Note that fiH/n'u = Ay /Ay. We assume 
that photo-dissociation equilibrium holds for this second 
box, too, i.e. 

x'coUUw Sdust{A'y) Sn^iNli^) Sco{Nco) ^ n'^^U^-jx', 

The question is now for which value of the extinction 
A'y (or equivalently for which mean density n'^j) do both 
boxes have the same CO abundance x'qq = xco? Once 
we know A'y such that x'qq = xqo (without necessarily 
knowing xco or x'qq) we can then use © to obtain xqq 
and, hence, x'qq. Since N'^q = x'^qu'^V cx A^x^q/Z' 
and iVco oc Ayxco/Z, we get iVco = [Av / A'y)N'^Q. 
The H2 column densities TVhj and are determined 
by the chosen H2 mean abundance model ^ and the 
values of Ay and A'y , respectively. 

In order to determine A'y s.t. x'qq = xco we demand 
that the solution reduces to A'y = Ay for Umw = C^mw 
and assume that sum involving the rate coefficients on 
the r.h.s. of pOl) and (fTTj) is not strongly affected by a 
change in the UV field. The latter assumption is based 
on a plausibility argument, namely that since xco = 
xco for A'y^ also the abundances of other species i likely 
satisfy x'^ ~ x^. 

Dividing ((TT]) by (HO]) we obtain 

1 ^ A'y S^^^Ay) SnANn,) ScojN'^^Ay/A'y) 
Umw Ay Sdust{A'y) Su2{Nii^) 5'co(-^co) 

(12) 

This equation has to be solved iteratively (see below) to 
obtain A'y and Nco = {Av / A'y)N'QQ for a giv en Ay and 
Umw- The shielding factors are tabulated in ILee et ahl 
(|1996D . see also appendix IB] 

CO self-shielding and H2 cross-shielding are due to a 
large number of individual absorption lines in the near 
UV. Hence, we multiply the column densities that enter 
the shielding functions by Lc/ L to take into account that 
the bulk velocity of the gas changes of the order of the 

^ Each shielding factor varies between and 1, where 5 = 
corresponds to perfect shielding. 
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Fig. 1. — Mass-weighted mean CO abundance xco a-s func- 
tion of mean extinction Ay. Symbols show the results of 
dr iven turbulence, magneto-h ydrodynamics simulations presented 
in IGlover fc Mac Low! I I2011I ). The star symbols correspond to 
a simulation without UV radiation, the other symbols are for 
'^MW = 1- Lines show the predictions of our mo del (see text) 
that is calibrated using the simulation results of Glove r Ik Mac Low! 
II2011I ). The lines correspond to simulations with different metal- 
licities: Z = Zq (solid black lines), 0.3 Zq (dashed magenta lines), 
and 0.1 Zq (dot-dashed red lines), and different UV radiation 
fields: Umw = 10"'', IQ-^, 1, 10^, 10* (set of Unes from top to 
bottom). Our model predicts that at fixed Ay the CO abundance 
varies strongly with the interstellar radiation field, but not with 
metallicity - except at high visual extinction, where the number 
density of carbon atoms in the gas, which scales with Z, imposes 
an upper limit on the CO abundance. 

velocity width of the shielding line over a finite coher- 
ence length Lc- We use Lc = 1 p c, similar to the co- 
heren ce length for H2 self-shielding (iGnedin fc Kravtsovl 
I2011D . to convert mean abundances into shielding col- 
umn densities. We show in appendix iBl that shielding of 
CO by dust grains is the dominant shielding mechanisms 
if Lc ^ 1 pc. Hence, none of our results depend on the 
precise value of Lc provided it is sufficiently small. 

If we ignore self-shielding and shielding by H2, and 
approximate the dust shielding factor as a plain expo- 
nential 5dust = e~'^°^^ , we arrive at a simpler version of 
equation ()12p that contains only the unknown Ay: 



1 



Umw 



(13) 



We first solve (fT^ with the exponent^ 70 = 3.4, and 
then use the obtained approximate solution as a trial 
value to solve (IT2t for A\/. 



In Fig. [T] we show the model predictions for the mean 
CO abundance for Umw = 10"^ - 10^- For a Milky 

® This value provides a reasonable approximation, for our pur- 
poses, of the dust shielding function over the range Ay = — 9 . 
We found that using a more sophisticated approximation in I I13I I. 
or using the tabulat ed s hielding function directly, does not affect 
the fina l so lution of I I12I I that we are trying to solve for. The solu- 
tion of JTSj approximates the final solution reasonably well, except 
at low metallicities Z < 0.1, since dust shielding is the dominant 
shielding factor for CO (see appendix [Bj . 



Way like radiation field {Umw — 1) the CO abun- 
dance changes strongly with Ay, but not with metallic- 
it y. The latter is by co n struct ion because the simulations 
of IGlover fc Mac Low! (|2011f ) do not show a significant 
metallicity dependence. A stronger (weaker) UV field 
lowers (raises) the abundance of CO at a given mean ex- 
tinction as expected. For Umw 1 the CO abundance 
depends (even at fixed Ay) on metallicity. This behavior 
is a consequence of H2 cross-shielding and can be under- 
stood from (I12p as follows. 

At a given Ay a larger UV flux lowers the CO abun- 
dance, hence Ay < Ay- The corresponding shielding 
columns Nn^ and A^^ are larger at lower metallicity 
(since Ay is kept fixed) and, due to the shape of the 
shielding function Sh^, the ratio S^^lN^^^) / Sii^iNii^) > 
1 is larger at lower metallicity. From (fT^ it can be seen 
that a larger ratio <S'h2(A^H2)/'^H2(^H2) can be compen- 
sated by increasing UmWi while Ay and Ay remain fixed. 
Hence, at a given Ay, lowering of the CO abundance re- 
quires a larger UV field in the low metallicity case than at 
high metallicity. Thus, raising the UV field {Umw > 1) 
leads to a larger CO abundance when the metallicity is 
lower. Analogous argument shows that lowering the UV 
field {Umw < 1) leads to a smaller CO abundance when 
metallicity is lower. If dust were the only shielding com- 
ponent, the CO abundance would only depend on Ay 
and not explicitly on metallicity. 

Combining ([SlSTf)) one gets (with T = 10 K): 



Tio ~ 5.9 X 10 



1.1 X 10^ 



-17 



cm 
3 km s~^ 



3 km 



Av 



Zr. 



z 



Av 

xco^v 



(14) 



The CO optical depth tiq increases rapidly with Ay due 
to the strong dependence of xqo on Ay. Specifically, 
for Av ^ 3K km s^^ and Umw ^ Ij '''10 varies between 
- IO-^Zq/Z, ~ Zq/Z, and ~ IO^Zq/Z for Ay = 1, 
3, and 10, respectively. We see that in this case the 
gas becomes optically thick to the CO line emission for 
Ay ^ 3, almost independent of metallicity. This is a 
significantly larger visual extinctio n compared w i th the 
results from the PDR calculation bv lWolfire et~all (poTof ) 
who predict that no = 1 should occur around Ay ~ 1. 

With Xco and x^.^ at hand, the column densities A^co 
and A^Hai and the optical depth of the CO emission line 
Tio can be derived ([5]). The CO velocity integrated in- 
tensity Wco can be calculated via ([2]) , if the brightness 
temperature and the CO line width is given. We discuss 
the brightness temperature in section 12.31 Our simula- 
tions do not resolve the velocity field within molecular 
clouds, hence we have to make a sensible choice for line 
width Av. We consider two possibilities. 

• The first case {constant line width) assumes that 
the velocity width is constant and has a value typ- 
ical of molecular clouds in the Milky Way Av = 3 
km s^^. This ansatz is motivated by the fact that 
our X-factor model is based on ISM simulations of 
a fixed box-size (20 pc) and 3 km s~^ is consistent 
with the observed size- Hnewidth relation of molecu- 
lar cloud s of such size (|Larsonlll981l : iSolomon et al.l 
[1987. : .Hever et al.l[2009ir 
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The second case [virial line width) assumes: 



Aw = 



-,1/2 



Fig. 2. — X-factor of the J = 1 — )■ ^^CO transition as function of mean extinction as predicted by our model (see text). (Left) Model 
predictions for a constant CO line width Av = 3 km s~^, (Right) for a virial scaling of the line width (see text). Symbols and lines are 
as in Fig.[T] Solid black lines show the predictions for different radiation fields {Umw = 10~^, 10~^, 1, 10^, Iff*) in increasing order from 
bottom to top. The horizontal, dashed line indicates the canonical galactic X-factor of Xqo.mw = 2 x 10 cm" -2 K-l km-l s. The dotted 
line indicates the trend of X-factor with Ay found in IGlov cr & Mac Low 11201 II ) for Umw = 1- The thin dot-dashed line in the bottom 
right shows the expected scaling Xco ■^v when the 1-0 line is optically thick and the line width fixed. 

density with Ay. In fact, the X-factor increases 
Unearly with Ay/Z, if we assume that the CO Une 
width and the excitation temperature remain fixed, 
and with the square- root of Ay jZ if we assume a 
virial scahng of the hne width. 

• The X-factor at fixed Ay increases with increasing 
UV field in the optically thin and moderately opti- 
cally thick regime of the CO line, e.g., for Ay < 4 
at Z = O.lZg and for Ay < 8 at Z = \Zq, but it 
is insensitive to the UV field in the highly optically 
thick regime. 

• For T = 10 K, Au = 3 km s'^ and a UV field in 
the range J7mw — 0.01 — 100 the minimum of the 
X-factor is at 0.7-1.1 x 10^" cm^^ Yr^ km'^ (for 
Z = Zq), 1.7 - 2.8 X 1020 cm-2 K^^ km^^ (for 
Z = O.SZq), and 4.5 - 6.3 x 10^" cm'^ R-^ km"! 
(for Z = Q.IZq). The minimum of the X-factor can 
be lowered further by increasing the metallicity to 
super-solar. 



which is based on assumption that massive molec- 
ular clouds are bound and cl ose to energy equipar- 
tition or even virialisation (iSolomon et al.l [19871 : 
lHeyer"eraI][200lt [Ballesteros-Paredesll2006D . Here, 
/grav is a constant of order unity that depends on 
the density structure and geometry of the cloud, 
its degree of virialization, and line-of-sight projec- 
tions. We assume /grav = 1- 

Equations ([1]), ([2]) imply that Xqo cx (for a con- 

1/2 

slant line width) and Xqq cx Sjj^ (for a virial line width) 
in the optically thick regime. These scalings bracket 
those computed from detailed rad iation transfer calcu - 
lations of MHD CMC simulations (jShettv et al.ll2011b[ ). 

In Fig.[2]we show the X-factor N^i^/Wco as function of 
mean extinction, for different metallicities and radiation 
fields, and for the case of a constant and virial line width 
scaling, respectively. The main predictions of the model 



• The X-factor depends on mean extinction Ay. It 
reaches a minimum at Ay = 2 — 10 (depending on 
metallicity and UV radiation field) , roughly where 
the CO line becomes optically thick (tiq « 1). It 
increases at higher and lower values of Ay. The 
increase at low Ay is a consequence of the much 
faster decline of the abundance of CO (compared 
to II2) with decreasing dust column. At large Ay, 
both hydrogen and carbon are fully molecular, the 
CO emission is saturated, and the X-factor in- 
creases simply due to the increase of the H2 column 



• The X-factor at fixed Ay does not strongly depend 
on metallicity at sufficiently low visual extinction. 
On the other hand, Xqq depends on metallicity at 
large Ay in our model because Wco is approxi- 
mately constant in the optically thick regime and 
a change in metallicity at fixed Ay corresponds to 
a change in N^i.^ (the gas is predominantly molec- 
ular), and thus in Xco- 

2.2. Simulations 

A detailed description of th e simulations can be found 
in lGnedin fc Kravtsov! ()201in . Here we repeat the main 
points for the convenience of the reader, following closely 
iFeldmann et al.. (,20 111) . 
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All simulations are run wi th the Eulerian hydrodynam - 
ics + N-body code ART (|Kravtsov et all [l997l f200l . 
that uses an adaptive mesh refinement (AMR) technique 
to achieve high spatial resolution in the regions of inter- 
ested (here: regions of high baryonic density). First, we 
ran an initial cosmological, hydrodynamical simulation 
down to z = 4. This simulation follows a Lagrangian re- 
gion that encloses five virial radii of a typical L^, galaxy 
(halo mass ^ 10"'^^ Mq at 2; = 0) within a box of 6 comov- 
ing Mpc//i. The mass of dark matter particles in the high 
resolution Lagrangian patch is 1.3 x 10® Af©. We adopt 
the following cosmological parameters r^mattcr — 0.3, 

VtK = 0.7, h = 0.7, ribaryon = 0.043, = 0.9. 

This initial, self-consistent simulation is consequently 
continued for additional ~ 600 Myr before it is analyzed, 
but now with metallicities and UV fields fixed to a spe- 
cific, spatially uniform value. At the end of each simula- 
tion the high resolution Lagrangian region contains one 
massive disk galaxy (halo mass ^ 4.2 x 10^^ Af©) and 
several smaller galaxies. The spatial resolution is ^ 60 
pc in physical coordinates. 

We have run a grid of simulations with six different 
metallicities Z/Zq =0.01, 0.03, 0.1, 0.3, 1.0, 3.0 and 
three different values of the interstellar radiation field 
C/mw = 0.1, 1, 100. We continued one of our simulations 
{Z/Zq = 1, C/mw = 1) for additional 400 Myr and found 
no significant changes in the ISM properties. This indi- 
cates that the predictions of our simulations should also 
hold for redshifts z < 3, at least unless/until ISM prop- 
erties change radically. In fact, in the computation of 
the brightness temperature of the CO line we use the 
z = CMB temperature in order to compare the model 
predictions with observations in the local universe. We 
have also run additional simulations at solar metallic- 
ity and C/mw — Ij but with 2 times better and 2 times 
worse peak spatial resolution confirming that none of our 
results suffer strongly from numerical resolution effects. 

The simulations include a sub-grid model for star for- 
mation, metal enrichment, supernova feedback (type 
la and type II) and cooling by metal lines and II2. 
The molecular hydrogen fraction is computed self- 
consistently, including a chemical network comprised of 
6 species and radiative tra nsfer of the UV continuum and 
the Lyman- Werner bands (jGnedin fc Kravtsovll20li[l . 

2.3. Modeling the X-factor on suhgrid scales 

We include the X-factor model in our numerical simu- 
lations as follows. 

First, we compute the mean visual extinction for a cell 
given its metallicity, density and an estimate of its line-of- 
sight depth, see ((T]). We then predict the mass weighted 
CO abundance xco based on Ay, metallicity and UV 
radiation field as described in ii2.1l The H2 mass frac- 
tion can either be taken directly from the cell (computed 
self-consistently with the non-equilibrium chemical net- 
work and radiative transfer within ART) or it c an be 
computed from the table of iGlover fc Mac Low! ()2011f ) 
(ignoring its dependence on UV field). While the first 
approach is our default method, we arrive at similar pre- 
dictions for the X-factor in either case, see § 13.11 Next, 
we convert the mean abundances xqq and into col- 
umn densities (see, eq. ^ by multiplying the abundances 
by the Sobolev-like length igob = p/|2Vp|. This pro- 
cedure leads to more reliable column density estimates 



compared with the use of the cell size as the line-of-sight 
depth and also avoids a few spurious artifacts related t o 
the numerical resolution, see IGnedin fc Kravtsovl (IMl . 
Finally, we compute the CO intensity and the X-factor 
of the cell via eqs. (lUlS]). 

The gas temperature T and the velocity width Ad, 
which both may vary in space and time, enter the predic- 
tion of the X-factor through equations ([2]), (g]), (0. The 
temperature that is measured in individual grid cells is an 
average temperature over the unresolved thermal struc- 
ture of the gas and thus not necessarily a good estima- 
tor of the interior temperature of a molecul ar cloud, for 
which observations indicate T ~ 10 K fe.g.. iHever et al.1 
I2009D . We therefore fix the gas temperature to T = 10 
K. Temperatures may increase in galaxies with strong 
star bursts and with densities high enough to allow for 
heat transfer between dust and gas via c ollisions (e.g. 
IStacev et al.l 119911 : INaravanan et al.l 120111 ). but explor- 
ing this possibility is beyond the scope of this paper. For 
At; we consider the cases of a constant line width and 
that of a virial line width scaling, see ^,2.1\ 

3. THE X-FACTOR ON THE SCALES OF GMCS 

3.1. Scaling with H2 column density 

Figure [3] shows results of applying the X-factor model 
described in the previous section to 12 numerical sim- 
ulations with constrained ISM properties ranging from 
Z = 0.03Zq to 3^0 and C/mw = 0.1 to 100. AU runs fol- 
low almost the same Xco ~ A^H2 relation at high column 
densities, but break off from this asymptotic relation at 
a specific H2 column density. The scaling of this asymp- 
totic relation is Xqq oh (constant line width) and 

1 /2 

Xco oc Sjj^ (virial scaling) as expected. Note that the 
dependence of Ay on metallicity via equation ([7|) is the 
reason why this figure looks very different from Fig. [51 

Figure [3] also shows that the U V radiation field has an 
significant impact on the X-factor. However, it should be 
pointed out that at fixed H2 column density the effect of 
a change of the UV field by a factor 1000 can be often be 
compensated by a less than factor 2 change in metallicity. 

Both the H2 column density at which the minimum of 
the X-factor occurs and the corresponding minimal Xco 
value decrease with increasing metallicity. More specifi- 
cally, for a simulation with Zq and a UV radiation field 
in the range C/mw = 0.1 — 100 Xco stays within a fac- 
tor of 2 of the galactic value Xco,mw = 2 x 10^° 
cm^^ km~^ s over a wide range of II2 surface densities: 
2.5 X 10^1 cm-2 < Nu^ < 2 x 10^^ cm^^ for a constant 
CO line width and 2.5 x lO^^ cm^^ < iVjj^ < 5 x 10^^ 
cm^^ for a virial scaling, respectively. Assuming the gas 
is fully molecular, the latter range in H2 column density 
corresponds to a range in total gas mass surface density 
(inch He) ~ 50 — 1000 Mq pc~^. This covers the range of 
surface densities of GMCs found in a variety of studies, 
including those that determine gas surface densities with- 
out the use of ^^CO emiss ion (ILa rson 1981; Hcvc r et al.l 
120091: IRoman-Duval et al.ll2010inHeiderman et al.ll2010l) 
and which are thus not biased due to variations in Xco ■ 

Hence, we conclude that the near constancy of the 
X-factor among the CMC population in our Galaxy is 
caused by the coincidence that for near solar metallic- 
ity A"co changes only weakly within the typical range 
of GMC surface densities. However, the figure clearly 
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Fig. 3. — X-factor vs H2 column density on ~ 60 p c scal es as a function of metallicity and UV flux. (Left column) model predictions 
using the tabulated H2 fractions of Glover & Mac Lowj I I2011I) . (Right column) using H2 fractions computed via the photo-chemical network 
within the simulation. (Top row) assuming a constant line width {Av = 3 km s~^), (bottom row) analogous predictions assuming a virial 
line width scaling (Ad oc S^/^)^ ggg jjgj] The curves correspond to 12 cosmological simulations with constrained ISM properties. The UV 
radiation field varies between Umw = 0.1 (dot-dashed lines), Umw = 1 (solid lines), and Umw = 100 (dashed lines). The metallicity varies 
(from bottom to top; see legend) between Z = 3Zq, Z = Zq, Z = 0.3Zq, Z = O.IZq, and Z = O.OSZq. Each curve connects the median 
Xco for a given H2 column density. The 16"" and 84th percentiles of the Xqq distribution for Umw = 1 is shown by the gray shaded 
areas. The galactic X-factor Xco,MW = 2 X 10^" cm~-^ km~^ is indicated by a solid horizontal line. Circles with error bars show 

results for the Perseus and Ophiuchus clouds in the Milky Way and should be compared with the Z ~ Zq lines IHeiderman et aLll2010l : 
their Table 5). 



shows that Xqo changes strongly with H2 surface den- 
sity if it is outside the range of ^ 50—1000 Mq pc^^, and 
for higher or lower metallicities even within that range. 
Generally, at any metallicity, the X-factor is expected to 
vary with H2 column density unless either the H2 sur- 
face density distribution in the ISM is very narrow (so 
that the dependence cannot be studied) or limited to the 
nearly flat part in t he Xgo — -/Vho curve (the Milky Way 
coincidence; also cf. iShettv et al.ll20libl ). 

While the dependence of Xco on metallicity is obser- 
vationally confirmed, see section [3T3l the increase of the 
X-factor with Nn^ at large column densities is not yet 
clearly established. Column density observations that 



are based on dust extinction maps should provide a good 
way to test these pred ictions. For instance, table 5 of 
IHeiderman et al.l ()2010( ) provides H2 column densities es- 
timates both from ^^CO and dust extinction maps. If we 
assume that the latter is a reliable measure of the true gas 
column density then the ratio of the two column densi- 
ties provides us with a measure of the X-factor in units of 
^co.MW- We include the X-factor derived in this way in 
FigureO A weighted linear regression of such Xco — 
correlatior0 results in an exponent of 0.53±0.14. Obser- 

^ We only include data with extinction based H2 column densi- 
ties in excess of 10^^'^ cm~^ in order to avoid including the "up- 
turn" of the Xco ~ relation at low column densities in the 
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vations thus indicate an increase of Xco with Nn^ , con- 
sistent with a CO hne width that scales as Av cx S^/^, 
although the data does not exclude a somewhat steeper 
scaling. 

Figure [3] also shows the predicted scatter of the X- 
factor at fixed H2 column density. This scatter arises 
due to the degeneracy between gas density and H2 frac- 
tion that result in the same H2 column density, but which 
lead to different predictions for xqq , see © . This scatter 
should be treated as a lower limit on the actual scatter 
on such scales. The scatter increases towards lower H2 
column densities and varies, for Milky- Way like ISM con- 
ditions, between 0.3 dex at A^H2 2 x 10^^ cm^^, 0.2 dex 
at 7Vh2 ~ 3 X 10^1 cm-2, and < 0.1 dex at N^.^ > 5 x 10^^ 
cm~^. This relatively small lower limit (< 0.3 dex) on 
the scatter of the X-factor is also consistent with the no- 
tion of a nearly constant X-factor among galactic GMCs. 

Finally, comparing the left with the right column of 
Fig. [3] we can check how the X-factor predictions depend 
on the chosen H2 modelin g. The left column use s the 
tabulated H2 fractions of Glover fc Mac Low! (j2011f) and 
ignores any changes with UV radiation field, while in 
the right column the H2 fractions and column densities 
are computed self-consistently with the photo-chemical 
network in ART. 

Differences in the H2 modeling lead to shifts along the 
45 degree line since both Xqo and Nb_2 depend on the H2 
fraction. Specifically, the use of AR T-based H2 fractions 
ys usin g the tabulated H2 fractions of lGlover fc Mac Low! 
(|2011f) results in shifts of ^ 0.2 dex. However, these 
differences are relatively minor and in any case do not 
lead to qualitative changes of the results in the paper. 

3.2. A parametrization of the X-factor 

The dependence of Xco on H2 column density as 
well as on Wqo can be captured reasonably accurately 
(to typically within < 30%, except at low metallicities 
and/or H2 columns) by simple parametrizations. In fact, 
an appropriate rescaling of Xco, -/Vh2 and Wco with 
metallicity and Umw removes most of the trends shown 
in Fig. [3] and results in a simple one-to-one relationship 
between the rescaled variables, see Fig. S) We approx- 
imate the relationships between the rescaled variables 
with a function of the form 

y = il + e-yx)°'il + e-^/^xf, 

which reduces to a power law with slope —a and /3 in the 
limit a; —7' and x — >■ 00, respectively. 

Specifically, in case of a constant CO line width, the 
scaling between Xco and N^^ is well approximated by 

y = (Z/Z0)O-8[/-^O5^co/6.O X lO^^cm-^ R-^ km"! s, 
X = (Z/ZQ)0-83t/-^"67VHj2.5 X lO^icm-2, 
a = 5.5, and /3 = 0.91. 

For a virial scaling of the CO line width we find that 

y={Z/Z^r^'U^^"''Xco/l-l X lO^W^K-^km-^s, 
X = {Z/Zq)°-^U^%°''N^J2.2 X lO^icm-^, 
a = 5.5, and /3 = 0.445 

fit, see Fig. [3] 



provides a good approximation. 

In a similar fashion, the scaling between Xco and Wco 
can be approximated by 

y = [ZlZof-'^'^Xco /6.0 x lO^^cm"^ j^-i j^j^^-i 
X ^ (Z/ZQ)°-°3VF'co/75Kkms"\ 
a = 0.78, and /3 = 11.5, 

when a constant CO line width is assumed. In case of a 
virial scaling we obtain 

y = (Z/Zo)° '*''Xco/l-25 X lO^^cm^^ K^^ km^^ s, 
X = (Z/Zo)''-'*VFco/30Kkms"\ 
a = 0.78, and /3 0.88. 

The parametrization of the relation between Xco and 
A^Ha makes it straightforward to include our Xco niodel 
in numerical simulations without the full modeling de- 
scribed in section 12.11 The relation between Aco and 
Wco may be used to convert observed CO intensities 
into H2 column densities. This requires that the obser- 
vations reach a spatial resolution of ~ 60 pc. 

These power law approximations do not provide the 
overall scaling of the X-factor with metallicity for an en- 
semble of molecular clouds with a variety of properties. 
For that we need to marginalize over the II2 column den- 
sity distribution (or the distribution of CO intensities) 
which itself may depend on metallicity. We will discuss 
this issue in more detail in the next section. 

3.3. Scaling with metallicity 

In Fig. Owe show how the X-factor scales with metallic- 
ity and UV radiation field on ^ 60 pc scales. A clear pre- 
diction of our model is that Xco increases with decreas- 
ing metallicity in the range Z ~ 0.1 — ZZq. We compare 
our predictions with observations o f the X-facto r base d 
on infrared (IR) dust emission by iLerov et"aLl ()2011l) . 
These observations seem to indicate a slightly steeper 
slope, but the deviation to our predictions is within 1-2 
standard deviations of the formal fit error and hence not 
statistically significant. Our predictions also agree well 
wi th the s l opes a nd n ormalizations f o und i n the studies 
bv lWilsonI ()1995[ ) and lArimoto et all (|1996[ ). 

Our model predicts that the X-factor decreases with 
increasing UV field at sub-solar metallicity. A large 
UV field suppresses molecular clouds with r elatively low 
II2 column densities ( )Feldmann et al.ll2011| ) and, hence, 
clouds with large X- factors. Note this happens despite 
the fact that the X-factor at fixed H2 column density in- 
creases with UV field (at least for Z > O.IZ©). This 
(moderate) UV dependence of the X-factor may con- 
tribute to the rat her low Xco value s in some of the 
GMCs observed bv lLerov et alTipoTTI ). 

Fig. [5] also shows estimates of the X-factor based 
on virial ma s ses fr om high resolution CO maps by 
iBolatto et al.l (|2008f ). These observations do not feature 
strong metallicity trends, possibly due to the fact that 
they focus on CO bright clumps and do not account fully 
for CO-dark molecular envelopes around those clumps. 
The scatter in these observations is very large. Inter- 
estingly, our X-factor model predicts a similarly large 
scatter. 

Our simulations therefore suggest that one should ex- 
pect significant variations in the X-factor even at fixed 
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Fig. 4. — Reseated X-factor vs reseated H2 column density (teft eotumn) and reseated CO intensity (rigtit eotumn) on ~ 60 pc seates 
as a function of metatticity and UV ftux. (Top row) model predictions of ttie X-factor assuming a constant line width (A?j = 3 tcm s~^), 
(bottom row) analogous predictions assuming a virial line width seating (Aii oc E^/^). The solid red tine shows an empirical approximation 
of the relation between the reseated quantities. It approaches power law behavior at tow and high H2 column densities and CO intensities, 
respectively, see text. The reseating of the X-factor, H2 column density and CO intensity with a metallicity dependent factor removes most 
of the explicit metallicity dependence. 



metallicity and UV field. However, it is important to 
point out that the scatter depends on the CO sensitivity 
limit, with higher sensitivity (i.e., a lower limit) leading 
to a larger scatter. This result can be easily understood 
from Fig. [3] or Fig. S] Lowering a sufhciently small sen- 
sitivity limit further will imply that more regions with 
lower Wco (and hence lower A^Ha and larger Xco) are 
included in the analysis, hence increasing the overall scat- 
ter. For instance, for Z = Zq, Umw = 1 and a Wco 
threshold of > 0.2 K km our simulations predict a 
scatter (defined as half the distance between the 16% and 
84% percentiles of log^o -'^^co) of 0.45-0.5 dex, while it is 
0.25-0.3 dex for Wco > 1 K km s'^ The scatter is not 
strongly dependent on the interstellar radiation field over 
most of the studied parameter range (C/mw = 0-1 ~ 100, 



Z/Z0 = 0.1-1). 

We can fit the increase of Xco with decreasing metal- 
licity with a pure power law, a dependence that is often 
assumed in the literature, 

logio ^co = ai logio(2'/Z0) + oq. (15) 

We provide the fit parameters in Table [TJ For instance, 
for Umw — 1, Wco > 0.2 K km s~^, and in case of 
a virial scaling of the CO line width the slope of the 
Xco — Z relation is -0.74. What is determining this 
slope? Clearly, the slope of the Xco — Z relation depends 
on the slope of the Xco — relation (see Fig [2|) and 
that of the Ay — Z relation (see below). 

In order to study the latter we show in Fig. [S] the 
(volume-weighted) probability distribution functions of 
the mean visual extinction Ay, the hydrogen column 
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Fig. 5. — Dependence of the X-factor on metallicity and UV radiation field on ~ 60pc scales. Tlie Xqo predictions are based on 
cosmological simulations with constrained ISM properties and assume a constant CO line width Av = 3 km s~^ (left panel) and a virial 
line width scaling Av oc S^/^ (right panel). The solid red, green and blue lines show the median X-factor of all cells above the CO sensitivity 



threshold, Wco = 0.2 K km s'^, 



for a UV radiation field of Umw 



0.1, Umw = li and Umw 



100, respectively. The light shaded 

areas in red, green and blue show the 16**" and 84*'^ percentiles of the X-factor distribution for Umw = 0.1, Umw = li a^nd Umw — ^00, 
respectively. Magen ta squares and green circles show X-factor measurements of individual molecular clouds by Lcrov et al. (20l J) and 
IBolatto et al.l 1)20081 ). respectively. Dot-dashed lines show observed scalings of Xqo with metallicity by Wilson (1995) and Arim oto et al.l 
II1996I). The CO i ntensity threshold Wco = 0.2 K km s^^ roughly corresponds to the 3 — cr intensity cut for the SMC in the sample of 
ILerov et al.l 1)20111 ). The intensity threshold is higher for other (more metal-enriched) galaxies in their sample, typically Wco ~ 1 K km 
s~^. Applying this higher threshold has little impact on the median X-factor for Z > 0.3^0, but it narrows the width of the X-factor 
distribution by ~ 0.2 dex. The horizontal line at A"co,MW = 2 X 10^" cm"'^ km~^ s corresponds to A"co,MWi the canonical value of 

the galactic X-factor. On ~ 60 pc scales there is no unique X-factor for a given Z and C/mWi but rather a broad distribution with a median 
that increases with decreasing metallicity and decreasing strength of the interstellar radiation field. Given the uncertainties and the scatter 
in the observational data, the predictions of our X-factor model are consistent with direct measurements of Xco in molecular clouds. 



TABLE 1 

The dependence of Xco on metallicity on scales of ~ 60 pc 



CO line width 


Wco limit 


Umw 


ai 


ao 


3 km s^^ 


0.2 


K 


km 


s'" 


-i 


1 


-0.87 


20.49 


3 km s^^ 


0.2 


K 


km 


s" 


-1 


100 


-0.79 


20.41 


3 km s^^ 


1 


K 


km 


s" 


-1 


1 


-0.81 


20.40 


3 km s~^ 


1 


K 


km 


s~ 


-1 


100 


-0.71 


20.39 


virial scaling 


0.2 


K 


km 


s~ 


-1 


1 


-0.74 


20.46 


virial scaling 


0.2 


K 


km 


s~ 


-1 


100 


-0.53 


20.40 


virial scaling 


1 


K 


km 


s~ 


-1 


1 


-0.56 


20.42 


virial scaling 


1 


K 


km 


s" 


-1 


100 


-0.43 


20.41 



Note. — The first three columns denote (1) the assumption about the scaling of the CO line width that enters our model, (2) the 
minimum CO velocity integrated intensity of a 60 pc scale resolution element in order not t o be excluded from the Xco distribution, and 
(3) the normalized strength of the interstellar radiation field. The parameters of equation l|15|l . i.e., first order fit parameters between Z 
and the median of the Xco distribution, are provided in the last two columns. The fit parameters are calculated using a least squares fit 
over the range O.IZq < Z < IZq. 



density A'h, and H2 column density Nb_2 of all ~ 60 pc 
resolution elements above the CO sensitivity limit 0.2 K 
km s~^. This figure demonstrates that (1) the median 
A'h increases with decreasing metallicity, (2) the median 
Ay decreases with decreasing metallicity, and (3) the 
peak in the H2 surface mass distribution coincides with 
the peak in the hydrogen surface mass distribution. This 
latter point is a statement of the fact that a large frac- 
tion of the gas that is detectable in CO is hydrogen in 
molecular form. However, it is noteworthy that there 
is a significant population of ISM regions with relatively 
low A^Ha /-^H (and hence low H2 mass fractions) that still 



make it above the CO sensitivity limit, especially under 
low Z and high C/mw conditions. 

A simple fit of the change of the median Ay with Z 
over the range Z/Zq =0.1 — 1 gives Ay oc z° '^^~°-^. 
In addition, a comparison of the median Ay in Fig. [5] 
(e.g.. Ay ^ 3 for Z — Zq, Umw — 1) with the Xqo - 
Ay relation. Fig. [21 shows that for such Ay the Xco — 
Ay relation has a negative slope. Therefore, when the 
metallicity decreases, the median Ay decreases and the 
median X-factor increases. 

A very different approach to the one presented i n this 
paper has been pursued bv lKrumholz et ahl (|2011l) . We 
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Fig. 6. — Probability distribution functions of the mean visual extinction Ay, the hydrogen column density Nn and the molecular 
hydrogen column density A'^h, on ~ 60 pc scales as measured in simulations of varying metallicity and radiation field (see legend). The 
solid black line shows the distribution of Nn (bottom axis) and Ay (top axis) of all ~ 60 pc regions with a CO velocity integrated intensity 
Wco greater than 0.2 K km s~^. The sharp fall-off at low column densities is a consequence of the Wqo sensitivity threshold. The dashed 
blue line shows the corresponding distribution of 2 X Nn^ (bottom axis). All column density distributions are normalized to an integral of 
unity over the plotted range. The vertical solid line indicates Nn = 6 X 10^^ cm~^ and the vertical dot-dashed line corresponds to Ay = 2. 
An arrow near the bottom axis shows the median of the Nn and Ay probability distribution. The median Ay decreases with metallicity 
roughly as cx: ^0-25-0.3 Qvgj. ^j^g considered metallicity range Z/Zq ~ 0.1 — 1. Consequently, the median column densities of molecular 
clouds increase with decreasing metallicity. 



call this ansatz for modeling Xqo, whic h is based on 
the re sults of PDR model calculations by iWolfire et alj 
(|2010( ). the K/W model and discuss it in more detail in 
the appendix. The X- factor in the K/W model is given 
by 

^co = Xco,oe^^^-/^^ (16) 

where Xco,o is a normalization constant (we use 10^" 
cm~^ km~^ s), AAy a (weak) function of metallicity 
of order unity, and Ay the mean visual extinction of the 
molecular cloud. An important property of this model is 
that it does not contain an explicit dependence on the in- 
terstellar radiation field and, furthermore, the predicted 
Xqo ~ -^v relation is almost independent of metallicity 
at fixed Ay- Of course, in order to de rive a Xgn — Z 
relation we need to specify Av{Z). In iKrumholz et aLl 
(|201lD the basic assumption is that all molecular clouds 
have very similar column densities, resulting in Ay (x Z. 
Under this assumption the K/W model predicts a very 



steep Xco — Z relation. If, however, we post-procesJ*^ 
our set of ART simulations with the ansatz (IT6l) , instead 
of using the X- factor model presented in section 12.11 we 
find a shallower dependence of the X-factor on metal- 
licity. Again, this can be understood based on the fact 
that our simulations predict a weaker than linear scaling 
of Ay with metallicity. In fact, if we use A'h ^ 6 x 10^^ 
cm~^ (for Z = Zq) and the scaling Ay cx i.e., 
values based directly on what we measure in the simu- 
lations, we find that ([T6|) results in a Xqo — Z relation 
similar to what we get when we post-process our ART 
simulation suite with the K/W model. 

We note that the metallicity dependence of the X- 
factor has nothing to do with the fact that there are 
fewer carbon atoms in a gas of lower metallicity. We 
demonstrated in Fig [T] that the CO abundance is a strong 

We convert the density and metallicity of a simulation cell 
into a mean visua l ext inction by multiplying with the Sobolev-like 
length Lsob, see i|2.3l AA\a is a function of metallicity only and 
computed as described in the appendix. 
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Dependence of the X-factor on mctallicity accord- 
K/W model (see appendix |A] Krumholz ct al. 2011, 
Many lines and symbols are as in Fig. [5] cx- 
The thick black lines (solid, dashed, dot- 
dashed, and dotted) show the predictions of the K/W model for 
different assumptions about the scaling of a clouds mean visual 
extinction Ay with metallicity. In particular, the black solid line 
assumes that hydrogen column densities of all molecular clouds 
(independent of metallicity) are A^h,o = 5.7 X lO^'^ cm"'^ (cor- 
responding to a nucleon column density AfH.o/0.76 = 7.5 X lO'^^ 
cm~'^). The dashed line assumes that clouds have hydrogen col- 
umn densities three times larger. The dotted line corresponds to 
the case that all molecular clouds have the same visual extinction 
Ay = 4, i.e., that A'^h oc Z~^. The dot-dashed lines shows Xqo 
based on a K/W model variant where solar metallicity molecular 
clouds have a column density of Afjj o ^^'^ the mean visual extinc- 
tion scales as Ay oc Z"'^*. Finally, the solid green line with crosses 
shows the result of plugging the K/W X-factor model into our set 
of ART simulations with Umw = 1 ^^'^ computing the X-factor in 
the same way as done for Fig. \5\ The outcome is a significantly 
shallower scaling of the X-factor with metallicity compared to mod- 
els that assume a constant column density of molecular clouds. In 
fact, it is reasonably close to the dot-dashed line over the range 
Z/Zq ~ 0.1 — 1 which uses the scaling of the median Ay with 
metallicity from the simulations, see Fig. [6] 

function of Ay, but (at fixed Ay) not an explicit function 
of metallicity. The only exception occurs at very large 
Ay, but then the abundance does not matter because 
the line is saturated. In our model the Xco ~ Ay rela- 
tion does (somewhat) depend on Z for large, fixed Ay 
due to the saturation of the CO intensity which implies 
Xco oc 7Vh2 oc Ay/Z. However, the Xco — Z rela- 
tion does not change significantly if this Z dependence is 
eliminated. Instead, as we have shown in this section, on 
GMC scales the dependence of the X-factor on metallic- 
ity is primarily a consequence of the metallicity scaling 
of the mean visual extinction of molecular clouds above 
a given CO detection limit. 

3.4. Implications for surface densities of GMCs 

Many galactic and extragalactic surveys assume a con- 
stant value of the X-factor, close to the galactic conver- 
sion factor Xco.MW = 2 x 10^*^ K^^ cm^^ km"-^ s, to 
predict II2 masses or column densities from -'^^CO data. 
This may introduce potential biases in the inferred II2 



column density distributions. 

To address this question we show in Fig. [8] the prob- 
ability distribution function of Nn^ on ~ 60 pc scales 
for simulations with different metallicities and UV field 
strengths. 

The solid black lines in each panel show the actual 
distribution of H2 column densities as measured in the 
simulations, i.e., using the H2 density within each 60 
pc resolution element and converting it into a column 
density by multiplying it with Lgob, see §2.31 The figure 
shows that an increase in the UV interstellar radiation 
field has a significant effect on the Nb_2 distribution. Low 
H2 columns are suppressed and the peak of the distribu- 
tion narrows and shifts toward higher H2 column densi- 
ties due to the environmental dependence of the HI to 
H2 transition on Umw- The hydrogen column density 
distribution is affected to a significantly smaller degree 
by a change in the radiation field. 

The dashed black lines in Fig. [5] correspond to the H2 
distribution after discarding all ~ 60 pc regions with a 
CO velocity integrated intensity below 0.2 K km s^^. 
This CO sensitivity cut removes preferentially lines of 
sight with intermediate and low H2 column densities, 
e.g., ^ 10^^ cm^^ for Z = 1, Umw = 1, and leads 
to a much narrower and peakier distribution. 

The blue dot-dashed and red dotted lines in Fig.[5]show 
the inferred H2 column densities = XcoM^Wco, 

i.e., those derived by multiplying the CO integrated 
intensity with the galactic conversion factor Xco,mw- 
Specifically, for a Milky- Way like ISM the inferred 
range[3 of H2 column densities 0.2 dex for a virial 
scaling, ~ 0.05 dex for a constant CO line width) is sig- 
nificantly smaller than the true width of the H2 column 
density distribution {'^ 0.7 dex). The peak position of 
the true H2 column density distribution is at ^ 1.3 x 10^^ 
cni^^, while the inferred distributions peak at ^ 1 x 10^^ 
cni^^. 

In contrast to the actual H2 column density distribu- 
tion, the inferred H2 column density distributions due 
not change if a CO intensity limit Wco > 0.2 K km 
is imposed, at least for > 4 x 10^^ cm^^. Hence, 

the much narrower range of the inferred H2 density dis- 
tributions and the bias towards higher column densities 
is not the consequence of such a limit. Instead, it arises 
due to the scaling of the X-factor with column density, 
as we now demonstrate. 

The inferred H2 column density is given as 



7>^obs 



Figure [3] shows that in the optically thick regime the X- 
factor can be approximated as 



X. 



CO 



X, 



CO.MW 



for a constant CO line width, and as 



X. 



CO 



CO,MW 



1/2 



10^-^cm 

for a virial scaling of the line width. Consequently the 



11 Measured as HWHM/^2 In (2), where HWHM is the distance 
from the peak to half the maximum at higher column densities. 
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Fig. 8. — Probability distribution function of the H2 column density on ~ 60 pc scales for simulations of varying metallicity and radiation 
field (see legend). The column density distributions are normalized to an integral of unity over the plotted range. In each panel the solid 
black curve shows the actual H2 column density distribution of all 60 pc regions in the simulation. The dashed black curve shows this 
distribution after discarding regions that have CO integrated intensities (according to our X-factor model with a virial scaling of the CO 
line width) below 0.2 K km s~^. The dot-dashed blue (virial scaling) and dotted red (constant CO line width) cur ves show the H2 column 
density distribution that is inferred from converting the CO emission as predicted by our X-factor model (see i|2.1l l into A'^h2 with the help 
of the canonical X-factor Xqq mw = 2 x 10^'' cm~^ km~^ s. We note that these inferred H2 density distributions due not change 

over the plotted range depending on whether resolution elements with Wqq < 0.2 K km s~^ are discarded, because any such element has 
^Yi2 ~ ''^CO.MW W'co < 4 X 10^^ cm~^. The increase of Xco with A'^H2 biases the inferred H2 column density compared to the actual H2 
column density if a constant MW-like X-factor is assumed. This results in an apparent narrowing of the H2 column density distribution 
and leads to a peak near lO'^^ cm~-^. 



inferred H2 column density is 



iV°bB ^ io22cm-^and N^, 



obs 



(10' 



22 2^ 



.1/2 



for a constant line width and a virial scaling of the line 
width, respectively. 

In the optical thin regime the X-factor raises steeply 
(well above Xco,mw) with decreasing column density 
and hence in general N^_^ < Nb_2- This effect is partic- 
ularly visible in high UV, low metallicity environments 
where the inferred (but not the actual!) H2 column den- 
sity distribution has a significant tail towards low H2 
column densities. 

The bias that we describe above adds another com- 
plication to the intense discussion of whether molecu- 
lar clouds have a "constant mean surface density" (e.g 



Kegell 119891: lScalo"1990 l: iVazguez-Semadeni et al 
ElmegreenT 200 2: Ballc steros-Paredes fc Mac Low 



1997 



2002 



Lombardi et al.li2010i : iHeiderman et al.ll2010l:lLada et al.l 
Observations, that are not based on ^^CO, and 



hence do not suffer from the mentioned bias, demon- 
strate that it is true in a weak sense, namely as the 
lack of a strong correlation of the mean surface den- 
sities of clouds with their sizes or masses, see Fig. 13 
Yet, a stronger and more controversial interpretation is 
that all clouds have very simi lar mean surface densities. 
In fact. iLombardi et al.l ()2010 ) find, based on extinction 
measurements of a small sample of molecular clouds, that 
the mean surface densities are constant with only 15% 
scatter if the cloud ar eas are defined by a fix ed extinc- 



tion threshold fb ut cf. IGutermuth et al 



20T1I) . Using a 



similar approach lleiderma n et al.l ()2010[ ) find that the 
scatter is ^ 30% if surface densities are measured above 
a fixed extinction threshold of Ay = 2. 

How do these observations fit together with the result 
shown in the top panel of Fig. [51 namely that of a rather 
broad distribution of H2 column densities, or with the 
observations of significant variations in observed GMC 
surface densities shown in Fig. ^ 
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Fig. 9. — Mean gas surface density of molecular clouds as a 
function of cl oud mass from observations using non-^^CO tracer. 
ILarsonl 1 119811 ) combines data from the literature, primarily ^^CO 
observations, available at the time. It is thus a rather hetero- 
geneous compilation and shows a large scatter, but no obvious 
trend of surface density with cloud mass. Rece nt ^'^CO sur- 
veys with larger statistics and better resolution (Heyer et al.l|2009l; 
[Roman-Duval et al. 2010) find a trend of increasing surface den- 
sity with cloud mass (or radius). These surveys also find significant 
variations in t he cloud surf ace density at fixed physical scale. The 
LTE masses of lHever et al.l lf2009) are multiplied by a factor two as 
suggested by the authors. Studies that derive gas surface density 
from extinction maps show reduced scatter if the surface densi- 
ties are measured above a fixed extinction threshold (Ay = 2 in 
IHeiderman et aLlEOlO t. 

A hint to a possib le solution is that the study of 
iLombardi et alj ()2010D finds an order of magnitude varia- 
tion in the enclosed mass (and consequently surface den- 
sity) if the cloud mass is measured within an aperture of 
fixed size (and not within a fixed Ay contour) . It is thus 
conceivable that the "constancy of the mean surface den- 
sity" is simply a matter of the identification methocQ of 
the molecular cloud and its characteristic properties. 

Hence, surveys that either do not properly resolve the 
clouds or measure masses at a fixed physical scale should 
find large scatter in the mean surface densities. Since we 
measure the surface densities on a fixed ^ 60 pc scale we 
indeed expect to see are rather broad distribution of H2 
surface densities. 

4. SPATIAL AVERAGING AND THE X-FACTOR ON 
GALACTIC SCALES 

Many extragalactic surveys use ^^CO observations to 
infer spatially averaged H2 column densities on ~ kpc 

Clouds in the MW have all rather similar metallicities 
and a fixed extinction threshold thus corresponds to a fixed 
threshold of the local surface density. Therefore, if molecu- 
lar clouds h ad an approximately self-similar (fractal) structure 
IjElmegreen fc Falgaronc 1 99^ or similar column density distribu- 
tions (cf., IBallesteros-Paredes et al.ll20ril) . then the mean surface 
density should scale with the local surface density of a given con- 
tour. This does not necessarily mean that the surface density of 
the entire region that is molecular is the same for each cloud. 



patches of galaxies or even for gal axies as a whole (e.g., 
iKennicuttillQQSl: iBigiel et al.|[2008l ). In most cases a sin- 
gle, constant conversion factor is assumed. In section 
^ 33.31 we demonstrated and discussed the dependence of 
the conversion factor on metallicity and interstellar ra- 
diation field on GMC scales. We now discuss these de- 
pendences on larger averaging scales. A related question 
that we want to address is by how much the X-factor 
can vary around its most typical value. Such variation 
arise from (1) the change of Xco with Nn^, and (2) the 
scatter of Xqo at fixed iVH2 due to the degeneracy with 
II2 fraction and hydrogen column density, see §3.11 

We compute the X-factor on large scales as the volume- 
weighted average^ of the small scale ~ 60 pc resolution 
elements, 



(X 



CO) 



{Wco) J 



dV 



i.e., the spatially averaged X-factor on scale I at a par- 
ticular point P in the simulation volume is computed as 
the ratio of the sum of the column densities and the CO 
integrated intensities of all 60 pc resolution elements 
within a box of extent I centered on P. 

Figure [10] shows the resulting Xqo as a function of 
II2 column density on 1 kpc and 4 kpc scales. Compared 
with Fig. [3] the spatial averaging on kpc scales and above 
reduces the variation of Xco with A^H2 i especially for 
low A'h2 • For instance, at solar metallicity and on kpc 
scales, the X-factor changes by less than a factor of two 
when A^H2 changes by 2 orders of magnitude between 10^^ 
cm~^ and 10^'^ cm~^. In contrast, on ^ 60 pc scales the 
corresponding change is a factor ~ 100. 

This demonstrates that the median (or the mean) of 
the Xco distribution becomes less dependent on Xh2 as 
one goes to larger scales, but it does not tell us much 
about the scatter of Xco on such scales. We therefore 
show in Fig.[Tl]the Xco distribution on both > kpc and 
on ^ 60 pc scales. We only include those 60 pc, 1 kpc, 
or 4 kpc ISM patches that have a CO velocity integrated 
intensity of at least 0.2 K km s^^. 

The figure highlights two important effects of spatial 
averaging. First, the Xco distribution on > kpc scales is 
considerably narrower than the one measured on scales of 
GMCs. Also the peak shifts toward higher Xco values. 
Secondly, the distribution of logj^g Xco becomes more 
symmetric and bears a closer resemblance to a normal 
distribution. Hence, even if Xco is not constant on small 
scales and can, in fact, vary over an order of magnitude or 
more, the spatial averaging ensures that variations of the 
X-factor on > kpc scales are much smaller. For instance, 
the X-factor varies by typically less than a factor ^ 2 
around its peak value on > kpc scales. 

In Fig. [12] we quantify the dependence of the X-factor 
on metallicity and UV radiation field on 4 kpc scales. 

The use of a volume average instead of an area average can 
be justified as follows. First, cloud self-covering is presumably 
relatively small as the size of each resolution element (^ 60 pc) 
constitutes a significant fraction of the scale height of the gas disk. 
Second, even clouds that do spatially overlap, e.g., in lines of sight 
edge-on through a disk galaxy, likely have a large enough velocity 
difference so that their CO intensities can be added. Since the 
former problem only arises in (near) edge-on views the results in 
this section may be safely interpreted in any case as predictions for 
sufficiently inclined (closer to face-on) views on disk galaxies. 
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Fig. 10. — X-factor vs H2 column density on galactic scales. (Top row) model predictions of the X-factor assuming a constant 
line width {Av = 3 km s~^). (Bottom row) analogous predictions assuming a virial line width scaling {Av oc T,^^'^), see i]2.1l 
(Left column) spatial averaging on scales of 1 kpc. (Right column) spatial averaging on scales of 4 kpc. The curves correspond 
to cosmological simulations with constrained ISM properties. The UV radiation field varies between Umw = 0.1 (dot-dashed 
lines), Umw = 1 (solid lines), and Umw ~ 100 (dashed lines). The metallicity varies (from bottom to top; see legend) between 
Z = 3Zq, Z = Zq, Z = 0.3^0, and Z = O.IZq. Each curve connects the median Xco for a given H2 column density. The 16"' 
and 84th percentiles of the Xco distribution for Umw ~ 1 are shown by light gray shaded areas (for Z = 3Zq and Z = 0.1^©). 
The standard deviation of the average X-factor is shown by the dark shaded regions. Both, percentiles and standard deviation, 
are only trustworthy if the given H2 column density bin contains 5 or more data points (indicated by filled circles). The galactic 
X-factor Xco,MW = 2 x cm ^ km ^ is shown by a horizontal solid line. The X-factor on > kpc scales depends 

primarily on metallicity and only weakly on the H2 column density or the strength of the UV interstellar radiation field. 

can be found in Table [21 Quantitatively, the power law 
exponent does not seem to be a strong function of spatial 
scale (compare with Fig.[T2]with Fig.|3]and Table[l]with 
Tabled. The value of the power law exponent depends 
on assumptions about the CO line width and the detec- 
tion threshold in CO intensity. We typically find power 
law indices ai for the median Xco ~ Z relation in the 
range - [-0.5,-0.8]. 

The power law slope can be understood from the re- 
sults presented in ^3.'2\ Since our galaxy models have 
a homogeneous metallicity distribution, the metallicity 
rescaling of Xqo, -^Ha and Wco that removes the most 



Similar to the results on CMC scales, see section [3T3l we 
find that ^co is strongly metallicity dependent. How- 
ever, the UV field (in the range Umw = 0-1 — 100) plays 
now no important role. This can be understood from the 
fact that a change in the H2 column density distribution 
hardly affects the X-factor, see Fig. [TUl The UV field 
becomes relevant only at at Z < O.IZq, due to the de- 
crease of the X-factor with increasing UV field at fixed 
H2 column density for a low metallicity ISM as discussed 
previously. 

The scaling of Xqo with metallicity is close to a power 
law. The parameters of a first order fit (in log-log space) 
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when measured on > kpc averaging scales compared with CMC scales. Also the peak of the distribution shifts toward larger values of Xqo- 
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Our X -factor model predicts a strong dependence of Xco metallicity, although less steep compared with the predictions of lUenzel et al'l 
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-0.52 (virial line width scaling), respectively. Fit parameters are provided in Table [2l 
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of the metallicity dependence on small scales also works 
on galactic scales. In fact, Fig. [10] shows that Xco is 
almost independent of A^H2 fixed metallicity. Hence, 
y ~ Z^'Xco (see ^3.2|) is approximately independent of 
both Z and Ajja and, therefore, ^co oc Z~^. This scal- 
ing argument predicts ai = —7 — 0.8 in case of a 
constant CO line width and oi ^ — 0.5 for a virial scal- 
ing. These expectations agree well with the ones ob- 
tained from a direct fit (Table [2]). 

The metallicity dependence that we find is shallower 
than the dependence expected from attributing offsets 
of z > 1 star forming galaxies in the ~ '^sfb. 
relation merely to the metallicity dependence of Xqo 
(jGenzel et al.|[20TTI ). which requires a power law index of 
~ [—1.3,-1.9]. This could imply that some part of the 
observed offsets is not due to Xqo variations, but, possi- 
bly, due to an actual deviation from the — '^sfr 
relation, as is obs erved more dramatically in mergers 
(IGenzel et al.ll2010l) o r in local galaxies with high specific 
star formation rates (|Saintonge et al.l [20111 ). A depen- 
dence of the i^cOa^i (and/or 1^003^2) to ^^COi^o 
conversion factor on metallicity may be an alternative 
possibility. 

5. SUMMARY AND CONCLUSIONS 

A proper physical interpretation of observational data 
based on CO emission requires a thorough understand- 
ing and modeling of the conversion factor, the X-factor, 
between molecular hydrogen and the J = 1 emis- 
sion line of carbon monoxide. We presented a novel ap- 
proach to study the properties of the X-factor based on 
a combination of high resolution (^ 0.1 pc) MHD simu- 
lations of the ISM (Glover & Mac Low 2011) as a "sub- 
grid" model and gas distribution of ^ 60 pc scales de- 
rived from self-consistent cosmological hydrodynamical 
simulations. This approach is feasible, because the spa- 
tial resolution of our cosmological simulations (~ 60 pc) 
is reasonably well matched to the box sizes of the ISM 
simulations 20 pc). In fact, this scale, which roughly 
corresponds to the scale of GMCs, is a natural scale to 
separate physical processes occurring on galactic and cos- 
mological scales from those relevant for star formation. 

The main advantage of our approach is that it takes 
into account the complicated interplay of chemical evo- 
lution, heating, and cooling in a turbulent ISM on small 
scales and is thus a step forward compared to, e.g., one- 
dimensional steady state PDR models or to ad hoc mod- 
els of small scale physics. A caveat is, of course, that 
it also suffers from any inaccuracies or missing physics 
in the underlying small scale simulations, but we are 
hopeful that these problems can be alleviated in the fu- 
ture. Our approach extents previous work that study 
the X-factor on ISM scales (e.g., iShettv et all 1201 lallbh 
to galactic scales and can be easily used in cosmological 
simulations, where all cell properties (densities, metallic- 
ities, radiation fields) are known, and thus provide self- 
consistent boundary conditions for the sub-grid model- 
ing. 

The main predictions of our model and the main results 
of this paper are as follows. 

1. The X-factor on CMC scales depends sensitively 
on metallicity, dust extinction, II2 column density, 
while two orders of magnitude variations of the UV 



radiation field lead only to moderate changes of 
Xco at a fixed hydrogen column. The changes in 
Xco are even smaller at a fixed H2 column density, 
because a higher UV field tends to reduce the II2 
column density at a given hydrogen column, which 
offsets to some extent the increase of Xco at fixed 
hydrogen column. 

2. In galaxies with solar metallicity the X-factor on 
CMC scales is predicted to be ~ 2 - 4 x lO^" 
cm~^ km~^ s (variations due to dependence of Xco 
on iVua), in agreement with the canonically as- 
sumed galactic conversion factor, for a wide range 
of gas column densities (^ 50 — 500 Mq) and with 
relatively small scatter at fixed Nn^ (< 0.3 dex). 
Our model also predicts an increase of the X-factor 
with decreasing metallicity that is in quantitative 
agreement with direct measurements of Xco from 
IR dust emission. 

3. The simulations predict an H2 column density dis- 
tributions with a peak around ^ 30— 100 Mq pc~^ 
for an ISM with solar metallicity, in agreement with 
GMCs observations in the Milky Way. 

4. We show that CMC column densities inferred from 
^^CO (J = 1 — >■ 0) observations via a constant 
galactic conversion factor are biased toward higher 
column densities and show a strong peak in their 
distribution, see Fig. \8\ In addition, such an ap- 
proach makes the spurious prediction of many lines 
of sight with low II2 columns in ISM regions that 
are exposed to a large UV field. These observa- 
tional biases are a consequence of relying on a con- 
stant galactic conversion factor that ignores the 
scaling of Xco with II2 column density. 

5. Spatial averaging, from CMC to ~ kpc scales (and 
beyond), decreases the dependence of Xco on the 
H2 (or total gas) column and the scatter in the 
X-factor. The column density dependence only re- 
mains at very low or large H2 columns, see. FigfTOl 
The UV radiation field plays only a small role on 
~ kpc scales, and hence, to a good approximation, 
the metallicity is the primary driver of the X-factor 
on > kpc scales. 

6. On kpc scales and above the X-factor scales with 
metallicity approximately as a power law with an 
exponent in the range [-0.5,-0.8], depending on 
assumptions about the CO line width and the ap- 
plied CO intensity detection threshold. The power 
law exponent on CMC scales is similar [—0.5, —0.9]. 

Recent CO emission studies (iDaddi et all l2010bllal: 
Genzel et all 120101: iTacconi et al.l 120101: lEmonts et al.l 
20 lit) have improved our knowledge of the relation be- 
tween star formation and its fuel, molecular gas, out to 
redshift 3. With instruments of the next generation, e.g., 
the square kilometer array, it will further be possible to 
detect CO bright gas at even higher redshifts. Going 
back in time these galaxies presumably differ strongly 
in their properties from their present-day (or even their 
z ~ 1 — 3) counterparts. Since all these studies observe 
CO we need to understand whether this gives unbiased 
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TABLE 2 

The dependence of Xco on metallicity on scales of 4 kpc 
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Note. — As Table [2] but for spatial scales of 4 kpc. 



estimates of the H2 density in the environments probed 
by the observations. Hence, a clear understanding and 
precise modeUng of the systematic trends of the X-factor 
with metahicity and H2 column density (among other 
factors) will remain a crucial challenge in order to 
properly interpret these future observations. 
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APPENDIX 
THE KRUMHOLZ / WOLFIRE MODEL 

In this section we give a short outline of the Krumholz/ Wolfire (K/W) X-fact or model a nd co mpare its predictions 
with those presented in section [2?T] The K/W model was introduced bv . Krumholz et al.l ()201li ) and is largely based 
on the results of PDR simulations by iWolfire et al.l ()2010( ) . 

The latter authors studied the case of a spherical cloud with an overall 1 /r density profile embedded in an isotropic 
radiation field. The cloud consists of cold (T < 300 K) gas clumps of density ric that contain most of the cloud mass 
and a warm (T ^ 8000 K) interclump medium with large volume filling factor. For a given mean visual extinction 
Ay, metallicity Z and impinging radiation field Umw Wolfire et al. compute (1) the radius Rco at which the optical 
depth of the J = 1 — > transition to the cloud surface is unity, and (2) the radius Rn^ at which the H2 mass fraction 
is 50%. Hence, they arrive at a prediction for the ratio between the cloud mass within Rn^ and the cloud mass within 
^co 

^^^^^^4 =exp(4 [0.53 - 0.045 ln(C/Mw/ric[cm-3])_ 0.097 ln(Z/ZQ)]/Av). (Al) 
IV1[< Rco) 

The depend ence of the mass ratio on Umw / ndcm^'^] can be eliminated with some additional modeling ([Krumholz et alj 
[2008l[2€09l) 

UMw/nc[cra-^] = 0.044(1 + 3.1[Z/Zo]°-365)/4.i. (A2) 
Based on the expression for the mass ratio it is now possible to construct a simple estimate for the X-factor 

^H. Mh. „ M (< i?H. ) M (< Rco) 



Xco = = — ^ = 0.76 , 

Wco Leo M{<Rco) L 



CO 



where Leo is the CO luminosity and the factor 0.76 accounts for the presence of Helium in the cloud. Observations 
of CO bright region s indicate that the l ast factor is roughly independent of the cloud environment, i.e. metallicity 
and radiation field ([Bolatto et al.|[2008[ ). but one should keep in mind that the observed scatter is large. Hence, 
for clouds that contain CO optically thick sub-regions (for optically thin clouds the last factor is undefined and the 
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Fig. 13. — X-factor of the J = 1 — ^^CO transition as function of mean extinction in the K/W model. The solid, dashed, and 
dot-dashed gray lines are the predictions of the K/W model for solar, 0.3 solar, and 0.1 solar metallicity, respectively. Filled circles and 
diamonds show the predicted X-factor for a nucleon column density A^n = Afjj/O.TG = 7.5 X 10^^ cm~^ (~ 60 Mq pc~^) and 1.5 X lO^'^ cm~^ 
(~ 120 M(7) pc~^), r espect ively. The other symbols show the results of driven turbulence, magneto-hydrodynamics simulations presented in 
IGlover Hi Mac Low! 1)20111 ) . The star symbols correspond to a simulation without U V radiation, the other symbols are for simulations with 
Umw = li see Fig. 12] The horizontal, dashed line indicates the galactic X-factor Jfco.MW = 2 X 10^*^ cm"'^ K"'^ km"'^ s. The dotted line 
indicates the trend of X-factor with Ay found in IGlover &: Mac Lowl (2011) for Umw = 1- The thin dot-dashed line in the bottom right 
shows the expected scaling Xqq oc Ay when the CO line is optically thick and the line width fixed. For a Milk y- Way like ISM (Umw = 1 ; 
Z = Zq) the predictions of the K/W model agree reasonably well with th e res ults of the ISM simulations bv IGlover fc Mac Low! I|20111 ) 
and, hence, with the predictions of our X-factor model presented in section [2. II except possibly at Ay < 0.4 and Ay > 10. However, the 
K/W model, which is not explicitly dependent on the strength of the UV field, differs noticeably fro m the predictions of our X-factor model 
for Umw 1- Also, it differs from the results of the no UV field run bv IGlover fc Mac Lowl 120111 ). 



method breaks) the X-factor should scale with M(< i?H2)A^(< Rco)- We fix the normalization by enforcing that 



Xco ^ -'^co,MW for a nucleon column density = Nii/0.76 — 10^^ cm ^ (~ 80 Mq pc ^) and Z — 



The K/W model, defined by equations (|A1[) . (|A2[) . and (jA3p . provides an estimate of Xqq as function of Ay and 
metallicity. 

Fig. [T51 shows the predictions of the K/W model. The X-factor depends primarily on visual extinction, i.e., at a 
fixed visual extinction the dependence on m etallicity is relatively weak . Quantitatively, at [/mw = 1, it agrees nicely 
with the results from the ISM simulations bv IGlover fc Mac Lowl ([20ITI ) over the visual extinction range 0.4 < Ay < 2 
and thus with the results of our model presented in Fig. [5] over such a range. A noticeably difference between the K/W 
model and our X-factor model appears at high visual extinction. Here, the K/W model predicts roughly constant 
X-factor at Ay > 3, asymptotically approaching ~ 10^° cm^^ K^^ kni^^ s, while our model predicts that the X-factor 
should depend on metallicity (the precise scaling depends on assumptions about the CO line width) . We note that the 
simula tions bv IGlover fc Mac Low! (I20T1I) are also indicative of such a metal-dependent scaling (see also iShettv et al.l 
I2011b( ). The other main difference between the K/W model and our model is that the X-factor in the K/W model 
is not explicitly dependent on the strength of the interstellar radiation field, while in our model the X-factor varies 
significantly with Umw at a fixed visual extinction, especially at Ay < 3. 

In order to provide Xco as a function of met allicity alone a further , and crucial, assumption about the visual 
extinction of molecular clouds has to be made. iKrumholz et alj ()2011|) assume (private communication) that the 
nucleon column density of all molecular clouds is ^ 7.5 x lO^'^ cm~^ and, therefore, that the visual extinction of 
molecular clouds scales with metallicity Ay cx Z. Fig. [T3] explains why this assumption leads to a very steep increase 
of the X-factor with decreasing metallicity. If the column densities of molecular clouds were, e.g., two times larger 
('^ 1.5 X 10^^ cm~^) a two times lower metallicity would be required to reach the same X-factor, resulting in a less 
steep Xco ^ Z relation over the observed metallicity range Z ^ 0.1 — IZq. Hence, the slop43 of the Xco — Z relation 
as predicted by the K/W model depends on assumptions about the average column densities of molecular clouds. 



Strictly speaking, the Xqq — Z dependence in the K/W model 



, , , , , , 1 -ji I I over a limited range in metallicities. 

does not obey a power law, but it can be fit with such a relation 
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Fig. 14. — CO shielding factors that enter the X-factor model, i|2.1l via equation I I12II for four different ISM conditions: (from top left 
to bottom right) Z = Zq k Umw = I, Z = Zq &c Umw = 100, Z = O.IZq & Umw = 1, and Z = O.IZ© & Umw = 100. Shown are the 
shielding factors S for dust shielding (dashed blue line), H2 cross shielding (dot-dashed red line), CO self-shielding (dotted green lines) and 
the total shielding (product of the former three). 5 = (or {I — S) / S = 00) corresponds to complete shielding from the UV interstellar 
radiation field, 5 = 1 (or (1 — S) / S = 0) is the no shielding case. The horizontal solid black line indicates S = 0.5. The vertical dashed 
(dot-dashed) line indicates the visual extinction that corresponds to a dust optical depth (optical depth in the CO J=l-0 line) of unity. 
The overall shielding of CO is dominated by dust shielding. 



Also, it is not entirely obvious why typical column densities of molecular clouds should not increase with decreasing 
metallicities in order to compensate to some extent for the loss in shielding. Observations of constant mass surface 
densities in CO bright clumps, such as found by B olatto et al.l (j200 8). do not rule out this possibility since at low 
metallicity most of the column density may stem from CO dark molecular gas or even from atomic hydrogen surrounding 
those clumps. Consequently, if the typical visual extinction of molecular clouds detectable in CO does not scale linearly 
with metallicity, the X-factor will increase more gradually with decreasing metallicity. For instance, in the extreme 
case that all CO detectable clouds have similar visual extinction, the K/W model would predict that Xco varies only 
very weakly with metallicity. 
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SHIELDING FUNCTIONS 

The X- factor has been measured by iGlover fc Mac Low! ()2011| ) for Umw — 1, but not for any other non- vanishing 
radiation field. The model presented in arrives at an estimate of the CO abundance for Umw 1 based on 
the assumption of photodissociation equilibrium of CO. This approach, equation (|12|) . takes various mechanisms into 
account that shield CO from the interstellar radiation field. Here we show these shielding factors explicitly for different 
ISM conditions and as a function of the mean visual extinction Ay ■ 

We compute for each Ay, Z, and Umw the CO a bundances using the model in and the H2 abundances according 
to the fitting formula in lGnedin fc Kravtsovl ()2011l ). Gas densities are obtained from Ay and Z assuming a line-of-sight 
length of 20 pc. We then convert GO and H2 abundances into shielding column de nsities by multip lying them with 
the coherence length Lc = 1 pc. Finally, we use the tabulated shielding functions of iLee etaP (|1996l ) to calculate the 
corresponding shielding factor as function of Ay. Note that the CO self-shielding factor is only known a posteriori, 
i.e., after equation is solved for the GO abundance. 

Fig. [T4l shows the different shielding contributions. For a sufficiently small coherence length {Lc < 1 pc), the shielding 
is dominated by dust shielding. 
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